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Abstract 

The magnetic susceptibility and Edwards-Anderson order parameter q of the spin-glass-like 
(SGL) phase of the double-exchange model are evaluated in the weak-coupling or RKKY limit. 
Dynamical mean-field theory is used to show that q = M(T/Tsgl)^i where M is the classical 
Brillouin function and T^gl is the SGL transition temperature. The correlation length of the SGL 
phase is determined by a correlation parameter Q that maximizes Tsgl and minimizes the free 
energy. The magnetic susceptibility has a cusp at Tsgl and reaches a nonzero value as T — > 0. 

PACS numbers: 75.40. Cx, 75.47.Gk, 75.30.-m 
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One of the most important models of itinerant systems, the double-exchange fDE) model 
■s believed to describe many physical systems ranging from the manganites Q to dilute 

magnetic semiconductors In the weak-coupling limit, the DE model becomes equivalent 
to a RKKY model with competing antiferromagnetic (AF) and ferromagnetic (FM) 
Heisenberg interactions between classical spins at every site. Recent work Q| has revealed 
that these competing interactions can stabilize a phase with short-range but not long-range 
magnetic order. We now show that the magnetic susceptibility of this spin-glass-like (SGL) 
phase has a cusp at Tsgl, marking the onset of short-range order, and reaches a nonzero 
value as T 0. The Edwards-Anderson (EA) order parameter g p,|a] of the SGL phase is 
identical to the square of the classical Brillouin function M(T/Tsgl)- 

The DE model contains a kinetic term that describes the hopping of electrons between 
neighboring sites and a potential term that aligns the electronic spins with the classical local 
moments at every site. There is no quenched disorder in the DE Hamiltonian 

H = -tY, {4aCja + c]aC»a) " 2 J2 " ^i, (l) 

where c|q, and Cjq, are the creation and destruction operators for an electron with spin a 
at site i, Sj = {l/2)cl^aai3C^ij is the electronic spin, and Sj = Snii is the classical spin of 
the local moment. Repeated spin indices are summed. When Jh > 0, the Hund's coupling 
favors the alignment of the local moments with the electronic spins. Due to the electron-hole 
symmetry of the DE model, we shall only consider electron concentrations p between and 
1 carriers per site. For small p, the hopping of electrons between neighboring sites favors the 
alignment of the local moments and the FM phase is stable. Since electrons with parallel 
spins cannot hop between singly-occupied sites due to the Pauli exclusion principle, the AF 
phase is favored over the FM phase near p = 1. But as shown in Ref.jJ], the competing FM 
and AF interactions may actually favor a SGL phase over the ordered phases for small Jr. 

Developed in the late 1980's by Miiller-Hartmann [3] and Metzner and Vollhardt jq|, 
dynamical mean-field theory (DMFT) exploits the momentum independence of the elec- 
tronic self-energy in infinite dimensions. Even in three dimensions, DMFT is believed to 
capture the physics of correlated systems including the narrowing of electronic bands and 
the Mott-Hubbard transition j^. Within DMFT, the local effective action on any site is 
parameterized by a Green's function that regulates the hopping of correlated electrons from 
other sites. We will use DMFT to study the DE model on a Bethe lattice in infinite di- 



mensions. The bare density-of-states of a Bethe lattice with 2; S> 1 nearest neighbors is 



No{n) = {A/TTW)yl - {2fi/Wy, where W = A^^t is the bandwidth and /i is the chem- 
ical potential. Since the Bethe lattice is not translationally invariant, q = is the only 
well-defined wavevector [J] 

In infinite dimensions, the high-temperature non-magnetic (NM) phases of the Heisenberg 
and DE models have a vanishing correlation length ^. The SGL phase is a bulk solution of 
the DE model with some of the same characteristics as conventional spin glasses: a finite 
local magnetization and spin-spin correlations that decay exponentially over distance jlol |. 
The SGL phase is characterized by a correlation parameter Q, defined as the average over all 
neighbors of sin^(^^j/2), where di is the angle between the central spin and a neighboring spin. 
Overall, the neighboring spins describe a cone with angle 2 arcsin(-\/Q) around the central 
spin. The FM and AF phases have, respectively, Q = and 1. The magnetization about 
every site decays exponentially with a correlation length ^ = —a/ log \ '2Q — 1|, where a is the 
lattice constant. Notice that ^/a diverges in the FM and AF phases but vanishes in the NM 
state obtained by setting = 1/2. Mathematically, the SGL phase was first introduced by 
Chattopadhyay et al. pj|, although its physical significance was not recognized until later. 
In lower dimensions, the SGL phase evolves into the phase with incommensurate correlations 
obtained in Monte-Carlo simulations jl^ . 

The transition temperature TsglIp, Q) of the SGL phase may be evaluated from coupled 
Green's function equations !4]. In the weak-coupling limit J^S <^ W and T ~ {Jy^SY /W , 
TsGhip, Q) is implicitly given by the expression Tsgl(p, Q) = Jes{p, Q)/3 where 

JMp,Q) - -2iJuSn2Q - DTE ^ _ (2) 

with Rn = (-Zn + \Jzl - 1^74 ) /2, Zn = + /i, and Vn = (2n + 1)ttT. Since T <^W, the 
sum TJ2nF{i^n) is equivalent to the integral (l/27r) / dvF{v) and Jcs{p,Q) is independent 
of temperature. The relation Tsgl(p, 0} = Jesip, Q)/3 correctly reduces to the FM {Q = 0) 
result first derived in Refs. Q and 14 1 . A similar derivation of Tsgl{p,Q) was recently 
provided in Ref.|15j]. Of course, Tsgl{p,Q) vanishes in the NM state with Q = 1/2. 

After Tsgl(p, Q) is maximized with respect to Q, Tsgl(p) exceeds the Curie and Neel 
temperatures in the concentration range 0.26 < p < 1. The correlation parameter Q changes 
discontinuously at Tsgl(p) from 1/2 in the NM phase above to a value less than or greater 
than 1/2 in the SGL phase below. The ground state is NM with Q = 1/2 for a single 




FIG. 1: The Brillouin function M(r) and EA order parameters g(r) versus r = T/Tsgl- 

concentration close to p = 0.5. 

The temperature dependence of the local SGL order parameter M = |(mj)| on site 
i may be evaluated in a local environment fixed by the correlation parameter Q by 
integrating the local action over the Fermion variables. The probability for rrij to 
point at an angle cos 6' with respect to the local quantization axis is proportional to 
exp(MJeff/3cos 6*) = exp((3M cos 6')/r), where r = T/Tsgl- Consequently, M has the solu- 
tion M(r) = coth(3M/r) — r/(3M), which is just the Brillouin function in the S* oo or 
classical limit plotted in Fig.l. The result for the FM order parameter is not surprising 
considering the weak-coupling equivalence between the DE model and a Heisenberg model 
with RKKY interactions between classical spins. What is surprising is that the short-range 
order parameters of the SGL phase are identical to the long-range order parameters of the 
FM and AF phases. For small r, M(r) ^ 1 - r/3 - + ^^(r^). 

To zeroth order in J^S/W, p is given in terms of yU by p = l/2+{l /it) ^6y/l — 5|, 
where 6 = 2fi/W. Carefully accounting for the dependence of the chemical potential 
fi{p) on J-nS/W for a fixed p, we have generalized the T = relation derived in Ref. ^ 
for the energy difference AE{p,Q) between the SGL and NM phases: AE{p,Q)/N = 
— (3/2)A'f(r)^ Tsgl(p, Q), where both sides are evaluated analytically to order {JuSY /W . 
This relation is precisely the same as the MF result for the FM phase of a classical Heisen- 
berg model. By integrating the specific heat to obtain the entropy, we have formally con- 
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structed the free energy difference AF(p, Q) [17| between tlie SGL and NM pliases. Because 
{l/N)dAF/dQ = — (3/2)M(r)^ (9Tsgl/<9Q, the free energy is minimized by the same corre- 
lation parameter Q that maximizes the transition temperature. 

We now evaluate the magnetic susceptibility of the FM and SGL phases in the weak- 
coupling limit by utilizing the formalism developed by Fishman and Jarrell Treating 
the electronic susceptibility as a 2n^- dimensional matrix, the local-moment susceptibility as 
a scalar, and the cross terms as 2r;,^-dimensional vectors in Matsubara space, x(q = 0,iu!m) 
can be written as a (2n^ + 1) x (2n^ + 1) supermatrix. The total susceptibility x(q = 0, iujm) 



is obtained by taking — > oo and summing x(q = 0,iu!m] 
fixed external frequency Um = 2rmTT. As discussed in Ref. 
Bethe-Salpeter equation 



over all matrix elements for a 



18j . x(q = 0,iuJm) satisfies the 



X(q = 0,iuj„i) = = OJ^m) + X^°Hq = 0,iu;m)L{i^m)x{<l = 0,icj„), (3) 

where T{iuJm) is the vertex function and x^^K^ = 0, icUm) is the bare susceptibihty. Within 
DMFT, momentum conservation at the internal vertices of irreducible graphs is disregarded 
so that internal Green's functions are replaced by their local values. Consequently, T_{iujjn) 
is independent of momenta and may be evaluated from an identical Bethe-Salpeter equation 
where x(q = Oji^m) and x^^Hq = O,^^^^) are replaced by local susceptibilities at a site i. 
Because the total spin + Sj) is conserved, the q = susceptibility x(ct = 0,^0;^) is 

proportional to 6m,o- So we shall henceforth take Um = to evaluate the elastic susceptibility 

x(q = o). 

The magnetic field is taken to lie along the magnetization direction in the FM phase but 
is averaged over all orientations in the SGL phase, which has no net magnetization. For 
example, the bare local-moment susceptibility on site i is given by (3S'^^{mf^) — (mjQ,)^| =^ 
(3S\1 - g)/3, where g = in the SGL phase and q = 3M^ " 2(1 - T/Tc) > (T < Tc) 
or (T > Tc) in the FM phase. Following the derivation in Ref.[l8|, we find that the total, 
elastic susceptibility is given by 

Wq = 0) = gl +^ (4) 

Xiq uj 3 T-Tc(l-g) 2 ' ^ ' 

where Ses = 5* -|- JuSNQ^fi). Since No^fi) oc l/W, the electronic contribution Jf^ SNn ifi) is 



much less than the local-moment contribution S to Ses in the weak-coupling limit 



19|. The 



electronic contribution enlarges or diminishes Ses depending on the sign of Jh- The final 




FIG. 2: The elastic, q = susceptibility versus T/Tsgl for several electron concentrations p. 

term in Eq.(jD), Nq{^)/2, is just the electronic Pauli susceptibility. Since it does not depend 
on Jh and does not diverge as T — Tsgl, we shall neglect the Pauli susceptibility in the 



.(jlD with S^s = S IS 
I . Bare in mind that 



subsequent discussion. For a FM, the Curie- Wess susceptibility of 
precisely the same as the MF result for a classical Heisenberg model 
the broad analogies between the (D)MF theories of the DE and classical Heisenberg models 
only exist in the weak-coupling limit and disappear once Jh'S' becomes of order W . 

Using the low-temperature behavior of M(r), we find that gpM — > 1 — and x(q = 
0) ('5cfT/9)T/T^ as T in the FM phase. By contrast, gsGL 1 - 2r/3 and x(q = 
0) (S'^g/3)/(3TsGL/2 — Tq) as T ^ in the SGL phase. Because the SGL phase has no 
long-range order and the local moments have no preferred orientation for any Q between 
and 1, the zero-temperature susceptibihty does not vanish as p ^ 0.26 and Tsgl — *■ 
Tq. The magnetic susceptibility in the SGL phase is plotted versus r for several different 
concentrations in Fig. 2. As expected, the SGL susceptibility has a cusp at Tsgl, which 
develops into a divergence asp ^ 0.26 and Q — ^ 0. Notice that the normalized susceptibility 
Tsglx(ci = 0)/'S'eff vanishes as Q ^ 1/2 and Tsgl — * in the vicinity of p = 0.5. 

Comparing Eq.(0]) with the parameterization of Sherrington and Kirkpatrick (SK) j^, we 
conclude that q is the EA order parameter We have plotted gsGL and gpM versus r in 
Fig.l. Our result gscL = (m*)^ should be compared with the SK expression q = {{Si)'^)j, 
where the inner expectation value is a thermal average for a given set of exchange couplings 
and the outer is an ensemble average over the distribution of exchange couplings. For the 
DE model, the EA order parameter is the same at every site even without an average over 
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quenched disorder. The DMFT result for the DE model differs from the SK result for the 
random S = 1/2 Ising model in at least one important respect: the SK prediction for ^ is 
not identical to the 5=1/2 Brillouin function whereas ^/q^GL is equivalent to the S = oo 
Brillouin function. 

To summarize, we have evaluated the total, magnetic susceptibility and EA order pa- 
rameters of the SGL phase in the weak-coupling or RKKY limit. We find that the SGL 
susceptibility has a cusp at Tsgl and reaches a nonzero constant as T — > 0, as expected for a 
phase with short-range magnetic order. The SGL phase is characterized by the short-range 
order of (rnj) and by the nonzero value of (nij)^ at every site. These were the original SK 
criteria p] for the existence of a SG. Unlike the SG ground state of the random Ising model, 
the SGL phase of the DE model can be studied analytically in the absence of quenched 
disorder. 
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AC05-00OR22725 with Oak Ridge National Laboratory, managed by UT-Battelle, LLC. 
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